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(57) Abstract: The invention is directed to a digital signal processing method and system for digitally processing a discretized in-put 
signal (12) provided in form of a plurality of discrete input signal values at corresponding discrete input time instants which are 
distributed non-equidistantly in time. The method comprises a conversion of the input signal (12) to a discretized output signal (16) 
representing the input signal (10) with discrete output signal values at discrete output time instants which are distributed equidistantly 
in time, wherein the conversion is carried out by a linear interpola-tion (14). The weight factors (w(i,ri)) of the linear interpo-lation 
are derived from an impulse response function of a filter (18). 



WO 2004/042553 



PCT/EP2002/012409 



METHODS, SYSTEM AND COMPUTER PROGRAM 
FOR DIGITAL SIGNAL PROCESSING 

5 

FIELD OF THE INVENTION 

The invention relates in general to the field of digital sig- 
lo rial processing and more particular to methods, a system and a 
computer program for processing a discretized input signal 
provided in form of a plurality of discrete input signal val- 
ues at corresponding discrete input time instants which are 
distributed non-equidistantly in time. 

L5 

BACKGROUND OF THE INVENTION 

Digital signal processing (DSP) is the processing of signals 
20 by digital means. Historically the origins of signal process- 
ing are in electrical engineering, and a signal here means an 
electrical signal carried by a wire or telephone line, or by 
a radio wave. Following the development of signal processing, 
a signal nowadays is in general considered a stream of infor- 
ms mat ion, for example stock prices, data from a remote- sensing 
satellite, etc. (for a general introduction to signal proc- 
essing see the book: Introduction to Signal Processing, S. J. 
Orfanidis, Prentice Hall 1995) . 

30 In many cases, the signal is initially in the form of an ana- 
log electrical voltage or current, produced for example by a 
transducer. An analog signal must be converted into digital 
(i.e. numerical) form before DSP techniques can be applied. 
The conversion generates a discretized signal in the form of 

35 binary numbers. 
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In DSP discretized signals are commonly processed in a vari- 
ety of ways, for example by using a filter circuit in order 
to remove or at least reduce the unwanted part of the signal. 
DSP is often implemented in Digital Signal Processor chips - 
5 specialized microprocessors with architectures designed spe- 
cifically for the types of operations required in DSP. DSP 
technology is nowadays commonplace in such devices as mobile 
phones, multimedia computers, video recorders, CD players and 
automobiles . 

10 

DSP techniques often comprise Fast Fourier Transforms (FFT) , 
which allow the frequency spectrum of a signal to be rapidly 
computed. A variety of signal processing methods, software 
products and microprocessor devices applying such a FFT re- 

15 quire equidistantly-sampled input signals. Classical sampling 
techniques measure the amplitude of a continuous time signal 
at regular time intervals. In some special applications, the 
amplitudes of a continuous time signal are however not ob- 
tained equidistantly distributed in time. In this case, 

20 called event domain sampling, such standard techniques (e.g. 
FFT, etc.) are not applicable. 

In WO 01/87647 Al which shows the preamble of the claimed 
subject-matter, such an event domain sampling is disclosed in 

25 the context of measuring and analyzing the angular velocity 
of wheel axles in order to derive an automobile's tire pres- 
sure. To be able to apply classical signal processing tools, 
the document teaches to convert the event domain sampled sig- 
nal (in form of discrete input signal values at corresponding 

30 discrete input time instants which are distributed non- 
equidistantly in time) to a time domain sampled output signal 
(in form of discrete output signal values at corresponding 
discrete output time instants which are distributed equidis- 
tantly in time, i.e. with a constant sampling period). This 

35 conversion applies nearest neighbor linear interpolation. 
Furthermore, this document discloses a preprocessing of the 
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measuring signal in order to remove signal noise originating 
from sensor errors . 



5 OBJECT OF THE INVENTION 

The object of the invention is to improve the conversion of 
event domain sampled input signals to time domain sampled 
output signal . 

o 

SUMMARY OF THE INVENTION 

According to a first aspect the invention achieves the object 
5 by a digital signal processing method for digitally process- 
ing a discretized input signal provided in form of a plural- 
ity of discrete input signal values at corresponding discrete 
input time instants which are distributed non-equidistantly 
in time. The method comprises a conversion of the input sig- 
o nal to a discretized output signal representing the input 
signal with discrete output signal values at discrete output 
time instants which are distributed equidistantly in time, 
wherein the conversion is carried out by a linear interpola- 
tion. The weight factors of the linear interpolation are de- 
5 rived from an impulse response function of a filter. 

According to a second aspect the invention achieves the ob- 
ject by a digital signal processing method for digitally 
processing a discretized input signal provided in form of a 

3 plurality of discrete input signal values at corresponding 
discrete input time instants which are distributed non- 
equidistantly in time. The method comprises a conversion of 
the input signal to a discretized output signal representing 
the input signal with discrete output signal values at dis- 

5 crete output time instants which are distributed equidis- 
tantly in time, wherein the conversion is carried out by a 
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nearest neighbor linear interpolation. This interpolation is 
performed with a sampling frequency which is p times higher 
than the desired output sampling frequency, wherein p is a 
real number and the output sampling frequency is the recipro- 
5 cal of the constant differences between the discrete output 
time instants of the output signal. Subsequently, a digital 
lowpass filter with a cut-off frequency equal to half of the 
desired output sampling frequency is applied to the output of 
the nearest neighbor linear interpolation. Finally, the fil- 
10 ter output is decimated by said factor p. 

According to a third aspect the invention achieves the object 
by a digital signal processing system for digitally process- 
ing a discretized input signal provided in form of a plural- 

15 ity of discrete input signal values at corresponding discrete 
input time instants which are distributed non-equidistantly 
in time. The system comprises a conversion component for con- 
verting the input signal to a discretized output signal rep- 
resenting the input signal with discrete output signal values 

20 at discrete output time instants which are distributed equi- 
distantly in time, wherein the conversion is carried out by a 
linear interpolation. The weight factors of the linear inter- 
■ polation are derived from an impulse response function of a 
filter. 

25 

According to a fourth aspect the invention achieves the ob- 
ject by a computer program product including program code for 
carrying out a digital signal processing method, when exe- 
cuted on a computer system, for digitally processing a dis- 

30 cretized input signal provided in form of a plurality of dis- 
crete input signal values at corresponding discrete input 
time instants which are distributed non-equidistantly in 
time. The method comprises a conversion of the input signal 
to a discretized output signal representing the input signal 

35 with discrete output signal values at discrete output time 
instants which are distributed equidistantly in time, wherein 
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the conversion is carried out by a linear interpolation. The 
weight factors of the linear interpolation are derived from 
an impulse response function of a filter. 

5 Other features inherent in the digital signal processing 
method and system and the corresponding computer program 
product are disclosed or will become apparent to those 
skilled in the art from the following detailed description of 
embodiments and its accompanying drawings. 



DESCRIPTION OF THE DRAWINGS 



In the accompanying drawings , 

15 

Figs. la,b show two plots of a time domain sampled sig- 

nal (fig. la) and an event domain sampled 
signal (fig. lb) ; 

Fig. 2 shows a schematic block diagram illustrating 

20 the general method steps as well as the func- 

tional means of a preferred embodiment of the 
invention; 

Fig. 3 depicts the Epanechnikov kernel of the linear 

interpolation method; 
25 Fig. 4 schematically shows a rotational speed sensor 

for measuring the angular velocity of rotat- 
ing axles; 

Fig. 5 shows a flow diagram- of the algorithm for 

generating of the rotational speed sensor 
3 0 output; 

Fig. 6 schematically shows a realistic rotational 

speed sensor with sensor errors; 
Fig. 7 depicts a plot of the rotational velocity of 

an ideal toothed wheel (dotted line) and of 
35 an realistic toothed wheel with sensor errors 

(solid line) ; 
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Fig. 8 



shows the Fourier decomposition of an event 
domain sampled data set corresponding to the 
measurement of the angular velocity by a non- 
ideal rotational speed sensor; 
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Fig. 9 



Fig. 10 



shows a comparison of the Fourier decomposi- 
tions of a wheel speed sensor's output with 
and without a preferred signal preprocessing; 
shows a schematic diagram of the preferred 
signal preprocessing for removing harmonics 
which overlap with the signal of interest. 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 

15 For a better comprehension of the invention, some of the ba- 
sic terms used for its description are initially explained in 
more detail with relation to the sampling of a time signal. 
The invention is however not restricted to such sampled time 
signals, but may also use any other kind of signal which does 

20 not originate from a particular sampling procedure of a time 
continuous signal but has the same characteristics as de- 
scribed below. 

The classical sampling technique measures the amplitude of a 
25 continuous time signal y(t) at regular time intervals T s 



This series y k is herein referred to as time domain sampled 
30 signal. 

In event domain sampling the continuous time signal is sam- 
pled every time the amplitude of the signal passes certain 
predefined signal values y k 



y k =y(kT s ), 



fc = l,2,..„ 



N . 



(1) 



35 
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y(t k ) = y k , 



/c = l,2,... 5 



N . 



(2) 



thus generating a series of data pairs comprising a time in- 
stant t k and a predefined signal value y k . Besides, equation 
5 (2) implies problems when the continuous time signal is not 
monotonously increasing. However, in many practical applica- 
tions the continuous time signal is monotonously increasing 
and the predefined signal values are further uniformly dis- 
tributed, so that equation (2) can be formulated as 



where Y s is a constant amplitude sample period . 

15 Figs. la,b show a continuous time signal 1 sampled in the 
time domain 2 and the same signal sampled in the event domain 
3, respectively. In fig. la the time instants 4 at which 
measurements of the signal are performed, are equidistantly 
distributed over the time axes 5 with constant sample period 

20 t s . On the signal axes 9 are plotted the signal values y k ac- 
cording to equation (1) . Fig. lb shows the same time continu- 
ous signal, but event domain sampled whereby the predefined 
signal values y k are uniformly distributed with sampling pe- 
riod Y s according to equation (3) . It is clearly visible that 

25 the thus obtained time intervals are non-equidistantly dis- 
tributed on the time axes 8. 

The constant distance between two time instants in the time 
domain is called sample distance. The reciprocal of the sam- 

30 pie distance is called sampling frequency. Half of the sample 
frequency is defined as the Nyquist frequency and defines the 
maximum frequency obtainable in the discretized time domain 
signal (for an introduction to Fourier decomposition of time 
domain to frequency domain signal see the book of S. J. Or- 

35 fanidis) . The maximum frequency within the discretized time 
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y(t k ) = kr s 



k = \,2,...,N , 



(3) 
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domain signal is just half of the Nyquist frequency. 



In the following, preferred embodiments of the invention will 
be explained in general with reference to fig. 2. 

5 

According to fig. 2 an input signal 10 is provided in form of 
an event domain sampled input data set 12 comprising a plu- 
rality of discrete input signal values x i at corresponding 
discrete input time instants which are distributed non- 
lo equidistantly in time. This input data set 12 is converted by 
an interpolating means 14 to a time domain output data set 16 
representing the input signal 10 with discrete output signal 
values y n at discrete output time instants which are distrib- 
uted equidistantly in time. 

L5 

The input signal 10 may originate from a time continuous sig- 
nal (for example a measurement signal or a control signal) 
whereby the term "time continuous" means that a signal value 
exists for every time instant, hence the signal value is a 

io continuous function of time. When such a time continuous sig- 
nal will be measured (sampled) by a measuring device, it will 
be provided in form of the input data set 12 . Alternatively, 
the measurement device may directly output discrete signal 
values (like the wheel speed sensor described below) . The in- 

is put data set 12 thus may either comprise a list of value 
pairs of discrete signal values together with corresponding 
time instant values or, as in the special case of equation 
(3) , a list of discrete time instants only. 

io The linear interpolation means 14 carry out a linear interpo- 
lation of the event domain sampled data 12 and output corre- 
sponding time domain sampled data set 16. It uses weight fac- 
tors for the linear interpolation which are derived from an 
impulse response function of a filter 18. 

$5 
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In a preferred embodiment, the linear interpolation is car- 
ried out according to the following formula: 

y n = *(cO = J<^^ ■ (4) 

/ i 

Here, T iut is the sample period of the output data set 16, 
denotes the input signal value at the time instant 
denotes the output signal value at the time in- 
stant t° ut =nT int and w(i,n) denote the weight factors. The indi- 
ces i and n denote the i-th element of the input data set 12 
respectively the 22-th element of the output data set 16. In 
general, the index i of the sum runs over all elements of the 
input data set 12. In a specific implementation, however, the 
range of the index i will be limited in order to reduce com- 
putational complexity. In some preferred embodiments, the 
weight factors are normalized according to ^w(i 9 n) = 1 . This 

i 

will guarantee an unbiased interpolation if x(tj tt ) is con- 
stant . 

Preferably, the weight factors are chosen according to the 
following formula so that the linear interpolation is a con- 
volution interpolation: 

(5) 

Here, the function h{t) is the impulse response function of a 
continuous time filter 18. In order to reduce the computa- 
tional complexity a fixed number M of normalized weights 
might be preferably used. The output data values are then ob- 
tained according to: 



WO 2004/042553 PCT/EP2002/012409 

y(nl\ nt jO = — — - wfr ) . ( 6 ) 

Limiting the number of weights to M has the consequence that 
only a fixed number of components of the discretized impulse 
5 response of the filter is used. Preferably, the underlying 
filter is causal, requiring that t k <nT int . The upper limit k 
of the summation is preferably chosen in such that nT Ult -t k is 
minimal and positive. Furthermore, M shall be chosen suffi- 
ciently large, in order to adequately realize the underlying 
io filter function. This is achieved by choosing M such that 

\h(nT ini -t k _ M )\<s, \fm>M, (7) 

is fulfilled for a given s which is chosen sufficiently 
15 small. 

In summary, the above three requirements for the summation 
result in the following procedure. The computation of an ac- 
tual output data value y{nT int ) starts with the multiplication 

20 of the input data value x(t k ) at the time instant t k <nT int ly- 
ing nearest to the output time instant nT htt with the corre- 
sponding weight factor w{k,n) which is computed according to 
equation (5) (showing the relation between the discrete 
weight factors and the continuous impulse response function 

25 h(t)) . Then, it is checked whether the condition 

(K w ^-'*)l <s (c P- equation (7)) is fulfilled for a given s. 
If this condition is not fulfilled, a further multiplication 
corresponding to the "subsequent" index k-1 is carried out 
and the result is added to the result of the multiplication 
30 corresponding to the "preceding" index k. These steps are re- 
peated until the above condition is fulfilled. If the condi- 
tion is fulfilled, the sum of all the preceding multiplica- 
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tions is outputted as the output data value y(nT int ) - 



Furthermore, in a preferred embodiment the impulse response 
function of a lowpass filter having a selectable cut-off fre- 
quency is used. The advantage of such an underlying filter 
model is that the cut-off frequency can be chosen arbitrarily 
to reduce undesirable aliasing effects which may occur in 
classical interpolation methods as described in WO 01/87647 
Al due to the following effect. Since event domain sampled 
signals can contain arbitrarily high frequencies, the conver- 
sion into time domain sampled data leads to a loss of fre- 
quencies when the sampling frequency of the time domain sam- 
pled data set is too low (i.e. the sample distance is too 
large) . Merely increasing the sampling frequency however is 
in practice no remedy since it increases computer memory us- 
age and computational complexity. Furthermore, all frequen- 
cies of the event domain sampled input data set above the 
sampling frequency are folded in the frequency range of the 
output data set thus leading to a deterioration of the output 
data set. This effect is called aliasing. Using a lowpass 
filter for the interpolation as described above advanta- 
geously removes the frequencies above half of the sampling 
frequency and thus suppresses the aliasing, if the cut-off 
frequency is substantially equal to half of the desired sam- 
pling frequency. 

One suitable lowpass filter is for example the second order 
Butterworth filter, which may be formulated as: 

H(s) = ~ — with a x =42 and a 2 = 1 , (8) 

wherein H(s) is the frequency function of the Butterworth 
filter in the frequency domain, which translates into the 
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following impulse response in the time domain: 



h(t) = ^-e *t sin(^0 , t > 0 , (9) 

fo V2 



5 in equations (8) and (9) , f Q is the cut-off frequency of the 
Butterworth filter . 

In another preferred embodiment, the weights of the linear 
interpolation are chosen according to an Epan e chn i ko v filter 
■0 kernel (i.e., impulse response function) (see also V.A. Epan- 
echnikov, Non-parametric estimation of a multivariate prob- 
ability density, Theory of Probability and Its Applications, 
14:153-158, 1969), which may be formulated as: 

"'W-Srt 1 -?), (10 » 

wherein d denotes the size of the neighborhood (see below) , 
and Q + denotes the positive part of the function. The Epan- 
echnikov kernel is depicted in fig. 3 for d = 1. The size d 
o of the neighborhood defines the number of weights used for 
the interpolation of the output signal value y(nT int ) , and 
might be defined via the following surrounding I n of the out- 
put time instant nT im : 

5 In = {i\nT int -d<t t <nT int +cl] (11) 

The discretized version of the Epanechnikov kernel may be 
written as: 
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The effect of this particular weight factor selection is a 
smoothing linear interpolation of the input data set attenu- 
ating or even removing large frequencies. The parameter d can 
be referred to as the bandwidth of the filter, and is related 
to the cut-off frequency of the filter. Therefore, a possible 
choice for the parameter d is : 



In order to reduce the computational complexity of the inter- 
polations as described above, the weight factors of the in- 
terpolation can be pre-calculated and stored in a look-up ta- 
ble. 

According to a further embodiment of the invention, the 
weight factors w{i 3 n) of equation (4) are chosen such that the 
linear interpolation is a nearest neighbor linear interpola- 
tion. In the nearest neighbor linear interpolation, the 

y n z=y(t° ut ) for t k _ x <nT int <t k are obtained via the formula: 



Here, y and 1— y are the weight factors for the nearest neigh- 
bor linear interpolation. Since the input data values x(t k ) 
are event domain sampled, the weight factors y and 1-y will 
vary for every output data value y{nT int ^j . 

When converting the input data set 12, some input data values 
of the input data set 12 may be not be considered in the 
nearest neighbor linear interpolation 14 when the sampling 
period of the output data set 16 is too large. The frequen- 
cies contained in the input data set 12 above the Nyquist 
frequency of the output data set 16 are aliased. In order to 
avoid this aliasing effect the nearest neighbor linear inter- 
polation step is performed with a sampling frequency which is 



int 



(13) 




int 



(14) 
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a factor p times higher than the desired output sampling fre- 
quency (i.e., l/T int ) . Hereby, the factor p should be chosen 
large enough so that the thus increased sampling frequency 
now comprises all frequencies contained in the input data set 
12. In a subsequent step, a digital lowpass filter with a 
cut-off frequency equal to 1/ (2T int ) is applied to the result 
of the interpolation step. The lowpass filter eliminates all 
frequencies above half of the desired output sampling fre- 
quency. In a last step, the obtained intermediate data set is 
decimated by the factor p by keeping only every p-th element 
of the intermediate data set and throwing away the rest. 
Thus, the frequencies of the input data set 12 above half of 
the desired output sampling frequency are removed without 
aliasing effects. 

In a preferred embodiment, the event domain sampled input 
data set 12 which is converted to the time domain sampled 
output data set 16 is obtained by measuring the angular mo- 
tion of a rotating axle with a rotational speed sensor which 
is a typical example where event domain sampling occurs. The 
wheel axles in vehicles, the cam shaft in engines and the mo- 
tor axles in robots are some examples for the application of 
these sensors . 

A typical example for such a rotational speed sensor is shown 
in Fig. 4. This sensor comprises in principal an inductive 
transducer 3 0 (for example an inductive or optical trans- 
ducer, a Hall sensor, etc.,) cooperating with the teeth 32 of 
a toothed wheel 34. Every time a tooth 32 passes the trans- 
ducer 30 it generates a trigger signal. The rotational speed 
is determined using the time elapsed between two adjacent 
trigger signals. Between two adjacent trigger signals of the 
transducer 30, the toothed wheel 34 rotates around the angle 
a^2yv/N , where N is the number of teeth of the wheel 34. 
The trigger signals are sent to a microprocessor unit con- 
sisting of an internal counter with clock period T clocK which 



14 



WO 2004/042553 PCT/EP2002/012409 

generates the event domain sampled input data set 12 accord- 
ing to the algorithm depicted in Fig. 5. 



In an initialization step 4 0 a clock counter clock is ini- 
tialized to the starting value clock = 0, and in a subsequent 
step 42 an event counter k is initialized to the starting 
value k = 1. After waiting until the clock period T c3 . ock has 
passed in a step 44, the clock counter count is incremented 
in step 46. Steps 44, 46 are repeated until the next trigger 
signal is received. When in step 48 a next trigger signal is 
detected, the algorithm proceeds to step 50 where the current 
clock value clock is multiplied with the timer's clock period 
Fciock and the result is stored in t*, the ic-th input signal 
value of the time domain sampled input data set . The event 
counter k is then incremented in step 52 and the algorithm 
returns to step 44 waiting for a next trigger signal. The 
time domain sampled input data set is outputted by this algo- 
rithm in the form (f,,f 2 ,f 3 ,...) indicating the time instants at 
which the teeth 32 pass the transducer 30. In the following, 
the term rotational speed sensor refers to the rotational 
speed sensor device shown in fig. 5 in combination with the 
microprocessor unit implementing the algorithm of fig. 6. 

The above rotational speed sensors realize an event domain 
sampling as described above which can be easily seen by re- 
placing Y s of equation (3) by the angular distance a. The 
time instants t k returned by the rotational speed sensor cor- 
respond to the angles a k =k-a which replaces y(t k ) in equation 
(3) . In the notation of the interpolation equation (4) this 
is equivalent to: tf=t k and a fax. Thus, in the time in- 

terval the rotational speed a> k of the axle to which 

the sensor is attached is given by 

In practice, the toothed wheel 34 of the rotational speed 
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sensor is non-perf ectly manufactured or may suffer wear dur- 
ing operation and thus shows deviations in the tooth widths 
and tooth distances. Fig. 6 illustrates a toothed wheel 34 
wherein the outline of an ideal wheel is drawn with a dotted 
line and the outline of a more realistic, non-ideal wheel is 
drawn with a full line. The ideal tooth distance is a and 8. 
is the -of f set error angle between tooth i and tooth i-1. Fig 
7 shows, assuming that the wheel 34 is rotated with constant 
angular speed co, the angular speed measured by an ideal 
toothed wheel sensor, which is constant in time (dotted 
line) , and by a real, imperfectly toothed wheel sensor, which 
is periodically fluctuating (solid line) with a periodicity 
N. N is the number of teeth of the rotating wheel 32 (N = 5 
in the example of fig. 7) . 

Fig. 8 shows the Fourier decomposition of the event domain 
sampled data set corresponding to the measurement of the an- 
gular velocity of a non-ideal rotational speed sensor. The 
periodic disturbance causes harmonics in the frequency domain 
with basic frequency equal to the revolution f requencyco 0 /(2;r) 
of the wheel 34. In fig. 8, the angular velocity is equal to 
co 0 =90rad/s resulting in a basic frequency of 14.3 Hz. All the 
harmonics are multiples f m = ma) 0 /(2flr) of the basic frequency. 
The harmonics may overlap with the signal of interest, and 
hence should preferably be removed. 

In the following, a preferred embodiment for removing the 
harmonics which overlap with the signal of interest is de- 
scribed with reference to fig. 10. These harmonics are re- 
moved in a preprocessing step before the event domain sampled 
data set 12 of the wheel speed sensor is converted into the 
time domain. If for example the signal lies in the interval 
[3 0 Hz, 60 Hz] , all harmonics falling in this particular in- 
terval shall be removed. For example, if the velocity is 70 
rad/s, the harmonics m = 3 , 4 and 5 should be removed. 



16 



WO 2004/042553 



PCT/EP2002/012409 



The preprocessing adopted for rejection of the harmonics is 
based on computed angular velocity ® k = a / (t k —t k _ x ) (calculated 
in step 52) between two consecutive time instants t k of the 
event domain sampled data set which is obtained from the ro- 
tational speed sensor 50. The co A between two tooth passing 
events are influenced from the rotational speed of the axle, 
from the sensor errors and from vibrations of the wheel. By 
exploiting the periodicity of the sensor errors, the a> fc can 
be expanded using a Fourier series expansion (provided in 
step 54) expressing these sensor errors: 

°>k Yu m a »i cos(2*r/cm / N) -h b m sm(2n km / N) (15) 

where co° = 27t/(t k+N/2 - t h _ Nn ) is the hypothetical angular velocity 
between two tooth passing events corrected for sensor errors 
(exploiting that 2 s /=°>- The Fourier coefficients a m and b m 
correspond to the harmonics f m . In order to estimate the 
Fourier coefficients a m and b m a fitting procedure, prefera- 
bly a Least Mean Squares method 5 6 may be used. With the es- 
timated values a m and b m it is then possible to subtract in 
step 58 the periodic sensor errors from the signal a> k : 

& k = fi>* - X„ a »> cos ( 2 ^^ / *0 + K sin(27r tan/N). (16) 

In order to reduce the computational effort, it is possible 
to restrict the preprocessing to the Fourier coefficients m 
of interest, e.g. the harmonics which overlap with the signal 
of interest. This is achieved by restricting the sum in equa- 
tions (15) and (16) to the values m of interest. 

Fig. 9 compares the Fourier decomposition of a wheel speed 
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sensor 1 s output in the frequency range from 0 Hz to 150 Hz 
with and without the above preprocessing. The dotted line 
represents the original signal calculated from the measured 

values a> k . The solid line represents (& k , the sensor error 
5 corrected signal with suppression of the sensor error induced 
harmonics m = 3 , 4 and 5 . As can be easily seen from this 
comparison, only the specified frequencies for these harmon- 
ics are removed from the signal. Naturally, the particular 
preprocessing of sensor signals may be applied irrespective 
10 of the subsequent specific interpolation as described above. 

Preferably, the above described method for converting event 
domain sampled data sets into time domain sampled data sets 
is applied to signals provided by rotational speed sensors 

15 attached to an automobile wheel axle which measure the angu- 
lar velocities of a tire. The tire velocity information is 
used for a determination of the automobile's tire pressures. 
This is accomplished by analyzing the vibration modes in the 
tire's rotational velocity. In general, the tire can be mod- 

20 eled by a damped spring with a resonance frequency. The reso- 
nance frequency of a tire typically depends on the tire's air 
pressure. This allows the tire pressure determination by 
measuring vibration resonances in the tire's rotational ve- 
locity. 

25 

Wheel speed sensors as described above are used for example 
in anti-lock braking systems (ABS) which provide a tire-road 
slip control and are a frequently integrated safety feature 
of an automobile. Thus, the rotational wheel speed sensors of 
3 0 the ant i- lock braking system can be used simultaneously for 
tire pressure monitoring and tire-road slip control. 

Further details of the tire pressure determination are shown 
in WO 01/87647 Al , the content of which is incorporated 
35 herein by reference. 
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According to an exemplified embodiment of the invention, the 
digital signal processing system may be any machine capable 
of executing a sequence of instructions that specify actions 
5 to be taken by that machine for causing the machine to per- 
form any one of the methodologies discussed above. The ma- 
chine may be an application specified integrated circuit 
(ASIC) including a processor and a memory unit. The instruc- 
tions may reside, completely or at least partially, within 
10 the memory and/or within the processor. 

In particular, the digital signal processing system may be 
implemented in the form of a computer system within which a 
sequence of instructions may be executed. The computer system 

is may then further include a video display unit, an alpha- 
numeric input device (e.g. a keyboard), a cursor control de- 
vice (e.g. a mouse), a disk drive unit. The disk drive unit 
includes a machine -readable medium on which is stored the se- 
quence of instructions (i.e., a computer program or software) 

2 0 embodying any one, or all, of the methodologies described 
above . 

The computer program product may be a machine-readable medium 
which is capable of storing or encoding the sequence of in- 
25 structions for execution by the machine and that cause the 
machine to perform any one of the methodologies of the pres- 
ent invention. The machine -readable medium shall accordingly 
be taken to include, but not be limited to, solid-state memo- 
ries, optical and magnetic disks, and carrier wave signals. 

30 

In the above description, for simplification, it is only re- 
ferred to the method case. The system case can be easily de- 
rived from the latter one by replacing at the appropriate 
places in the description the expression 'step' for the 
35 method case by the expression 'component 1 for the system 
case . 
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Furthermore, all other publications and existing systems men- 
tioned in this specification are herein incorporated by ref- 
erence . 

5 

Although certain methods, systems and products constructed in 
accordance with the teachings of the invention have been de- 
scribed herein, the scope of coverage of this patent is not 
limited thereto. On the contrary, this patent covers all em- 
10 bodiments of the teachings of the invention fairly falling 
within the scope of the appended claims either literally or 
under the doctrine of equivalents. 
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What is claimed is: 

1. A digital signal processing method for digitally process- 
ing a discretized input signal (12) provided in form of a 
plurality of discrete input signal values at correspond- 
ing discrete input time instants which are distributed 
non-equidistantly in time, comprising: 

converting the input signal (12) to a discretized output 
signal (16) representing the input signal (10) with dis- 
crete output signal values at discrete output time in- 
stants which are distributed equidistantly in time, 
wherein the conversion is carried out by a linear inter- 
polation (14) ; 

characterized in that the weight factors (w(i 9 ?i)) of the 
linear interpolation are derived from an impulse response 
function of a filter (18) . 

2. The method according to claim 1, wherein the linear in- 
terpolation (14) is carried out according to the follow- 
ing formula: 



x\t" 1 ) denoting the discrete input signal value at the 

discrete input time instant f" , y{t° lit } denoting the dis- 
crete output signal value at the discrete output time in- 



3. The method according to claim 2, wherein the weight fac- 



*(c")=E^>M'f), 




st ant t, 



- ? tlt , and w{i 9 n) denoting the weight factors. 



tors w(i P n) 



are chosen according to the following formula: 



n 



out 
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5 4 . 
5. 

10 

6. 

15 7 . 
8 . 

20 

25 9 . 
10 . 

30 



wherein h is the impulse response function of a continu- 
ous time filter (18) . 

The method according to anyone of the preceding claims, 
wherein the filter is a causal filter. 

The method according to anyone of the preceding claims, 
wherein the filter is a lowpass filter having a se- 
lectable cut-off frequency. 

The method of claim 5, wherein the filter is a Butter- 
worth filter. 

The method according to claim 5 or 6, wherein the cut-off 
frequency is automatically selected in dependence of the 
sampling frequency of the output signal (16) . 

The method according to anyone of claims 3 to 7, wherein 
only those of the weight factors ( w{i,n) ) are used for the 
interpolation (14) whose corresponding value of the im- 
pulse response function is greater than a predetermined 
lower limit. 

The method of claim 2, wherein the impulse response of 
the filter (18) is an Epanechnikov kernel (fig. 3) . 

The method of claim 9, wherein the Epanechnikov ' kernel 
(fig. 3) is chosen according to the following formula: 



wherein d is a filter parameter which is chosen in de- 




(10) 



pendence of the sampling frequency, and (•)+ denotes the 
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11. A digital signal processing method for digitally process- 
ing a discretized input signal (12) provided in form of a 
plurality of discrete input signal values at correspond- 
ing discrete input time instants which are distributed 
non-equidistantly in time, comprising: 

converting the input signal (12) to a discretized output 
signal (16) representing the input signal (10) with dis- 
crete output signal values at discrete output time in- 
stants which are distributed equidistantly in time, 
wherein the conversion is carried out by a nearest neigh- 
bor linear interpolation (14) ; 
characterized in that 

the nearest neighbor linear interpolation is performed 
with a sampling frequency which is p times higher than 
the desired output sampling frequency, wherein p is a 
real number and the output sampling frequency is the re- 
ciprocal of the constant differences between the discrete 
output time instants of the output signal; 

a digital lowpass filter with a cut-off frequency equal 
to half of the desired output sampling frequency is ap- 
plied to the output of the nearest neighbor linear inter- 
polation; and 

the filter output is decimated by said factor p. 

12. The method of anyone of the preceding claims, wherein the 
input signal (10) is provided by a rotational speed sen- 
sor (3 0,32,34) which in particular measures the angular 
motion of an axle. 

13. The method of claim 12, wherein the input signal (10) is 
preprocessed in order to remove periodical occurring sen- 
sor errors (8.) originating from non-ideal rotational 
speed sensors (30,32,34), by: 

fitting (56) the coefficients {a m ,b m ) of a Fourier series 
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expansion (54) on the basis of the harmonics (fig. 8) 
corresponding to the periodicity of the periodical occur- 
ring sensor errors (8.) to the input signal (10); 
subtracting (58) the Fourier series expansion (54) with 
the thus obtained coefficients (a m ,b m ) from the input 
signal (10) . 

14. The method of claim 12 or 13, wherein the input signal 
(10) corresponds to the angular motion of at least one 
tire of an automobile and is used for a determination of 
the automobile r s tire pressure. 

15. The method of anyone of claims 12 to 14, wherein the in- 
put signal (10) is provided by the rotational speed sen- 
sor of an anti-lock braking system of the automobile. 

16. The method of claim 14 or 15, wherein the tire pressure 
is determined by a wheel vibration analysis. 

17. The method of anyone of the preceding claims, wherein the 
weight factors {w(i,n)) of the interpolation (14) are 
stored in a look-up table. 



18. A digital signal processing system for digitally process- 
ing a discretized input signal (12) provided in form of a 
plurality of discrete input signal values at correspond- 
ing discrete input time instants which are distributed 
non-equidistantly in time, comprising: 

a conversion component for converting the input signal 
(12) to a discretized output signal' (16) representing the 
input signal (10) with discrete output signal values at 
discrete output time instants which are distributed equi- 
distantly in time, wherein the conversion is carried out 
by a linear interpolation (14) ; 

characterized in that the weight factors (w(i 9 n)) of the 
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linear interpolation are derived from an impulse response 
function of a filter (18) . 



19. A computer program product including program code for 
carrying out a digital signal processing method, when 
executed on a computer system, for digitally processing a 
discretized input signal (12) provided in form of a plu- 
rality of discrete input signal values at corresponding 
discrete input time instants which are distributed non- 
equidistantly in time, comprising: 

converting the input signal (12) to a discretized output 
signal (16) representing the input signal (10) with dis- 
crete output signal values at discrete output time in- 
stants which are distributed equidistant ly in time, 
wherein the conversion is carried out by a linear inter- 
polation (14) ; 

characterized in that the weight factors (w(i,n)) of the 
linear interpolation are derived from an impulse response 
function of a filter (18) . 
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